Electronic and Magnetic Properties of RMnOs/AMnOa Heterostructures 
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The ground state properties of RMnOa /AMnOa (RMO / AMO) heterostructures (with R=La, Pr, 
a trivalent rare-earth cation, and A=Sr, Ca, a divalent alkaline cation) are studied using 
a two-orbital double-exchange model including the superexchange coupling and Jahn- Teller lattice 
distortions. To describe the charge transfer across the interface, the long-range Coulomb interaction 
is taken into account at the mean-field level, by self-consistently solving the Poisson's equation. The 
calculations are carried out numerically on finite clusters. We find that the state stabilized near the 
interface of the heterostructure is similar to the state of the bulk compound (R, A)MO at electronic 
density close to 0.5. For instance, a charge and orbitally ordered CE state is found at the interface if 
the corresponding bulk (R, A)MO material is a narrow-to-intermediate bandwidth manganite. But 
instead the interface regime accommodates an A-type antiferromagnetic state with a uniform — 
orbital order, if the bulk (R, A)MO corresponds to a wide bandwidth manganite. We argue that 
these results explain some of the properties of long-period (RMO)^^^ (AMO)^^ superlattices, such as 
(PrMn03)^/(CaMn03)j^ and (LaMnOa)^^^ (SrMnOa)^^. We also remark that the intermediate states 
in between the actual interface and the bulk-like regimes of the heterostructure are dependent on the 
bandwidth and the screening of the Coulomb interaction. In these regions of the heterostructures, 
states are found that do not have an analog in experimentally known bulk phase diagrams. These 
new states of the heterostructures provide a natural interpolation between magnetically-ordered 
states that are stable in the bulk at different electronic densities. 

PACS numbers: 73.20.-r, 73. 21. Ac, 75.47.Lx 



I. INTRODUCTION 

Modern fabrication technology allows for the growth 
of multilayer structures that are nearly perfect at the 
atomic level, namely with minimal roughness, employing 
a variety of transition metal oxides (TMO). Due to the si- 
multaneous participation of several degrees of freedom in 
TMO, it is expected that the artificial multilayer struc- 
tures made of these materials will exhibit much richer 
physics than in conventional semiconductor heterostruc- 
tures. Indeed, recent studies have revealed some fasci- 
nating phenomena, such as the reconstruction of spin, 
charge, and orbital orders at the interface. ^'^ Many inter- 
esting properties have also been theoretically predicted 
using a variety of many-body techniquesi^i^'^ 

Among the several ongoing efforts, there is a con- 
siderable interest in the analysis of RMnOa/AMnOa 
(RMO/ AMO) heterostructures, where R=La, Pr, ... 
refers to a trivalent rare-earth, and A=Sr, Ca, ... is a 
divalent alkaline element. ^'Mi^iiSiiiii^ii^iii At low tem- 
peratures, the bulk RMO is in an A-type antiferromag- 
netic (A-AFM) state, which is an insulator, whereas the 
bulk AMO is in a G-type antiferromagnetic (G-AFM) 
state, that is also an insulator. Upon doping, the alloy 
Ri_;iA2;Mn03 ((R, A)MO) exhibits a variety of states 
depending on the doping concentration x, which con- 
trols the charge density in the alloy. However, the het- 
erostructure RMO /AMO could potentially behave differ- 
ently from its parent bulk compounds. For instance, the 



transfer of charge through the interface caused by the 
different Fermi energies, and concomitant different elec- 
tronic density concentrations, of the superlattice compo- 
nents causes a distribution of charge that it is not ho- 
mogeneous along the growth direction. Hence, several 
states may exist in different regions of the heterostruc- 
ture. While far from the interface the behavior must be 
similar to the one in the bulk compounds, close to the 
interface it may occur that phases very different from 
those of the superlattice components may exist due to 
the charge leaking through the interfaces. For instance, 
in the short-period LaMnOs/SrMnOs (LMO/SMO) su- 
perlattices the regime close to the interface exhibits fer- 
romagnetic (FM) metallic behavior, which is different 
from either LaMnOs or SrMn03-iii£>iiii^ More inter- 
estingly, as the number of Sr layers exceeds a critical 
value, the metallic behavior gives way to an insulating 
one, displaying a metal-to-insulator transition (MIT). 
This suggests that the electronic reconstruction at the 
interface has a crucial effect on the physical properties of 
the heterostructure Theoretically, considerable progress 
has been made in describing the ferromagnetism induced 
by this electronic reconstruction j^i^i^^'^^i^^i^^ A recent 
study by the authorsi^ focused on the MIT in the short- 
period LMO/SMO superlattices. There, a FM metaUic 
state with a dominant — orbital order was found 
at the interface. The insulating behavior of the super- 
lattice was explained as induced by Anderson localiza- 
tion of this quasi two-dimensional (2D) FM state. Note 
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that limited by computing power, most theoretical efforts 
based on computer simulations concentrate on superlat- 
tices with a small number of Sr layers. This is enough to 
understand the FM in the short-period superlattice, but 
may not be appropriate to study the properties of longer- 
period structures. Another limitation is that most theo- 
retical studies focus on the LMO/SMO superlattice. The 
corresponding bulk compound (La, Sr)Mn03 (LSMO) is 
a wide bandwidth manganite. But not much is known 
about the interfacial state of the RMO/AMO superlat- 
tice if (R, A)MO is a narrow or intermediate bandwidth 
manganite. 

In this manuscript, we use the two-orbital double- 
exchange model for manganites, which has been suc- 
cessfully applied to the study of bulk Mn oxides, to 
the analysis of RMO/AMO heterostructures. These het- 
erostructures are assumed to be grown along the (0, 0, 1) 
direction (c axis). We here also assume that the length 
of the heterostructure is long enough that the states at 
the two ends of the heterostructure resemble their bulk 
counterparts, i.e., an A-AFM at the RMO side and a 
G-AFM at the AMO side. This assumption allows us 
to focus on the only interface present in the heterostruc- 
ture. By using a relaxation method introduced in the 
following section, we can obtain the magnetic and elec- 
tronic properties of the ground state of the heterostruc- 
ture. Different from previous efforts that concentrated on 
the FM tendency in short-period superlattices, we find 
that in the heterostructure considered in this paper the 
state stabilized at the interface depends on the band- 
width of the corresponding bulk compound (R, A)MO at 
electronic density close to 0.5. If (R, A)MO is a narrow- 
to-intermediate bandwidth manganite, a state resembling 
the well-known CE-state of bulk manganites is found to 
be stabilized at the interface. Hence our calculations sug- 
gest a CE-like interface state in the PrMnOa/CaMnOs 
heterostructurci^ But if (R, A)MO is a wide bandwidth 
manganite, an A-AFM state is found at the interface in- 
stead. This is consistent with recent experimental results 
on the (LM0)„/(SM0)2„ superlattices.-^^ 

In addition, it is important to remark that one of the 
main results of our study is the observation of states close 
to the interface that do not have an analog in experimen- 
tally known bulk phase diagrams. These states arise as 
interpolations between, e.g., the A-AFM and CE-states 
that dominate in the bulk and interfaces, respectively. 
For instance, canting of the spins in the CE zigzag chains 
creates a novel "canted CE state" . Also the relative spin 
angle between adjacent CE planes can be different from 
those observed in the bulk. And in some occasions, ar- 
rangements of spins and orbitals were identified that do 
not have a clear bulk analog in experiments. 

This paper is organized as follows. In Sec. II, we intro- 
duce the two-orbital model and the numerical method 
used to obtain the ground state properties of the het- 
erostructure. The results for (R, A)MO corresponding to 
narrow-to-intermediate bandwidth manganites are pre- 
sented in Sec. III. Results for (R, A)MO corresponding 



to wide-bandwidth manganites are provided in Sec. IV. 
Finally, in Sec. V a discussion of our results is given, 
followed by conclusions. 

II. MODEL AND NUMERICAL METHOD 

To investigate the physical properties of the 
RMO/AMO heterostructure described in Sec. I, 
here the two-orbital model for manganites^ will be 
applied. This model has been widely used before 
to study the properties of bulk manganites, and the 
following assumptions are also widely accepted: (1) 
the t2g electrons are considered as localized and are 
described as classical spins with magnitude S = 3/2; (2) 
the Jahn- Teller lattice distortions are also assumed to 
be classical; (3) the Hund coupling between the t2g and 
Cg electrons is assumed to be infinitely large so that the 
eg spin is always parallel to the localized t2g spin at the 
same site. Based on these assumptions, the Hamiltonian 
that will be applied to the manganite heterostructure 
reads 

^ = - E (C''a,4c,/3 + H.c.)+ J2 ^^s.-s, 

<ij> <i,j> 

+ ^(<?!>, + Wi)ni + A Y^iQiiUi -f QaiTf -I- Qa.rf ) 

i i 

+ \ Y.^'^qI + qI + qI)- (1) 

i 

The first term of Hamiltonian Eq.(l) denotes the two- 
orbital double-exchange hopping term, a and [3 run over 
the two Cg orbitals d-^2_y2 (orbital a) and (i3z2_i.2 (or- 
bital b) of a Mn ion. Cia (c^q.) annihilates (creates) an 
Cg electron in orbital a at site i with its spin parallel 
to the localized t2g spin at site Sj. r denotes the ex- 
change direction, giving t'^'' = f^"- = Stf = St'^J' = 3io/4, 
= t^^ = -tf = -4"" = Vita/A, tl" = tf ^ 

t^"- = and 4^ — ^0 (^o is set to be the energy 
unit). The hopping amplitude is affected by the factor 
ilij = cos(^)cos(^) + sin(^)sin(^)exp[-i(.^i - (pj)], 
where 9i and ipi are the angles of the t2g spins in spher- 
ical coordinates. Here, we will assume that the hopping 
amplitudes are the same for electrons on both sides of 
the heterostructure. This may not be realistic given the 
possible mismatch of lattice constants between the two 
different compounds. However, since the superexchange 
coupling is more sensitive to the change of lattice con- 
stant than the hopping amplitudes themselves, as a first 
approximation we assume that the hopping constants 
keep the same value on both sides while the superex- 
change couplings may become layer dependent. The sec- 
ond term is the standard superexchange interaction be- 
tween nearest-neighbor (NN) t2g spins. Here the t2g spin 
Si = (sin 9i cos ipi, sin Bishiipi, cos Oi) has been normal- 
ized to a unit vector (the actual S'=3/2 magnitude of 



3 



the spins is absorbed in the superexchange couphng). To 
consider the effect of possible distortions from a perfect 
cubic lattice, two couplings are used: jf^J = J^'^ if i and 
j are NN sites in the same layer (with same z coordinate) , 
and jf^J — J^^ if i and j are NN sites belonging to two 
adjacent layers. In the third term, 1/)^ corresponds to a 
site-dependent Coulomb potential that originates from 
the charge transfer through the interface, and it is de- 
termined via the Poisson equation as described below. 
Wi denotes the work function on either side, which is 
determined by the positions of the chemical potentials 
in the corresponding bulk materials. More details on 
(jii and Wi will be discussed later in this section, n is 
the eg charge density. The fourth term stands for the 
electron-phonon coupling. The Qs are lattice distortions 
for the Jahn- Teller modes {Q2 and Q3) and breathing 
mode (Qi). r = (t^, ,t^) is the orbital pseudospin op- 
erator. The last term is the elastic energy of the lattice 
distortions considered here. The extra factor 2 for the 
breathing mode suppresses this mode as compared with 
the Jahn- Teller modes that are the most active. 

For any given i2g-spin and lattice configuration the 
above Hamiltonian can be solved by numerically diag- 
onalizing the bilinear fermionic sector. The ground state 
is approached via a relaxation technique: the optimized 
configuration of the oxygen coordinates and the t2g spins 
is determined by minimizing the total energy of the 
Hamiltonian Eq.([T|). This method is first applied to a 
4x4x4 lattice with periodic boundary conditions to 
estimate properties of the bulk compound on each side 
of the hctcrostructure. In this publication, we will fo- 
cus on the interface between an A-AFM and a G-AFM. 
Hence, we will adopt a set of parameters J^^ and A that 
gives A-AFM and G-AFM ground states at the limits 
of Cg electron densities n = 1 and n = 0, respectively. 
The same set of parameters is then used to calculate the 
ground state of the heterostructure. The hctcrostructure 
is defined on a 4 x 4 x 8 lattice (see Fig. [ij , with peri- 
odic boundary conditions in the xy plane {ab plane) and 
open boundary conditions along the z direction (c axis). 
The initial spin configuration is set to be A-AFM on one 
side of the heterostructure (4 layers) with electron den- 
sity n = 1, and G-AFM on the other side (4 layers) with 
n = 0. Regarding the relative spin direction between the 
A-AFM and G-AFM states, their relative angles 6 and 
(j) were allowed to take 16 equally-spaced values in their 
respective ranges, i.e. [0,7r] for 9 and [0,2tt] for </>, thus 
giving 256 possibilities. For each of these 256 possibili- 
ties, and a fixed set of couplings in the Hamiltonian, an 
independent optimization of the classical variables was 
made, namely the energy was minimized by solving 



V^.(i/2b> =0 



(2) 



self-consistently using the Broyden's method)2S where 
-0i = {di,^Pi,Qii,Q2i,Q3i)- At each step we keep the 
two end layers to be dominated by the A-AFM and G- 
AFM states, respectively. Namely, the dominant layer 
wavevector for the Fourier transform of the spin-spin cor- 




FIG. 1: Schematic representation of the heterostructure stud- 
ied in this paper. It is defined on a 4 x 4 x 8 lattice. The 
bulk material on one side of the heterostructure is in the A- 
AFM phase, while the bulk material on the other side is in 
the G-AFM phase. 



relations must be (0,0) and (7r,7r), respectively, otherwise 
the configuration is discarded. The rest of the layers, of 
course, may have different magnetic orders. The opti- 
mized spins and oxygen coordinates of the ground state 
then corresponds to the configuration with the lowest 
overall energy after this long optimization process. 

The charge transfer through the interface is taken into 
account via the self-consistent solution of the Coulomb 
potential For the heterostructure discussed here, in- 
volving manganites only, we assume Wa = Wq and 
make them equal to zero for simplicity. This approxima- 
tion is in agreement with results of previous discussions.— 
Then, the charge transfer through the interface is fully 
driven by the charge density difference, i.e., the charge al- 
ways transfers from the A-AFM side to the G-AFM side. 
To properly describe the charge transfer, the long-range 
Coulomb interaction among the mobile eg electrons and 
the positively charged ionic background must be included 
into the Hamiltonian as 



Coul 



ato ^ 



1 n,n 



riin 



2 \n 



(3) 



where a — e'^ /eato is a dimensionless screening parame- 
ter. For manganites, it is known that to is of the order of 
0.5 eV and the lattice constant a ~ 4A, hence a depends 
on the choice for the dielectric constant e. However, e is 
both temperature and frequency dependent and for this 
reason an accurate estimation of a is not well known. 
In this paper, the Hamiltonian Eq.Q is studied over a 
broad range of e values, 2 ^ e ^ 20, corresponding to 

0. 2 ^ a ^ 2. ri is the position vector of the Mn site 

1. n.i is the local electronic density at site i. stands 
for the effective positive charge density on the i-th Mn 
site arising from the background ions. Note that to sim- 
plify the model, we have already assumed that all the 
charges from the background ions are located on the Mn 
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sites. Therefore, for the RMO/AMO heterostructure, nf 
is fixed to 1 at the A-AFM side, and to at the other 
side to enforce the charge neutrahty. 

Equation [3] is solved at the mean-field level by intro- 
ducing the Coulomb potential 



(4) 



We then find that the Coulomb interaction in Eq.([3]) re- 
covers the third term in Eq.([T]). In practice, the Coulomb 
potential is determined at each step of the relaxation 
procedure by self-consistently solving the Poisson's equa- 
tion 



:{{r 



(5) 



Numerically the following discretization is applied: 
d'^4)/dx^ = 4)i+i - 24)i + d'^^jdy^ = ^i+y - 2(j}, + 
(j)i^y, and d'^cj)/dz'^ = 0^+21 - 20i+£ + cj)i. The open 
boundary condition is applied along the z direction (out- 
of-plane) but periodic boundary conditions are applied 
along the x and y directions (in- plane) . 

In this work, up to 4,000 iterations per set of couplings 
(an iteration here is defined as an update of the entire 
set of classical variables in the cluster) are used to obtain 
the optimized spin and oxygen lattice configurations of 
the ground state. The most typical number of iterations 
is approximately 1,000. At each iteration of this relax- 
ation procedure, up to 2,000 additional iterations at a 
fixed set of classical variables are used to solve the Pois- 
son's equation. Note that the diagonalization of fermions 
must be performed at each step in solving this Poisson's 
equation. Then, in order to find the optimized configura- 
tion for the ground state on a 4 x 4 x 8 lattice, typically 
approximately 10^ times the exact diagonalization of the 
256 X 256 matrix is necessary. This is very CPU time 
demanding. It is for these practical reasons that only the 
4x4x8 lattice is used here to study the properties of 
the heterostructure. 



III. RESULTS FOR NARROW TO 
INTERMEDIATE BANDWIDTH MANGANITES 

In this section, results for the RMO/AMO heterostruc- 
ture will be discussed, where (R,A)MO corresponds to 
a narrow-to-intermediate bandwidth manganite, such as 
(La, Ca)Mn03 and (Pr, Ca)Mn03. Therefore, the results 
presented here are expected to best describe the prop- 
erties of LaMnOa/CaMnOa or PrMnOs/CaMnOs het- 
erostructures. 



A. Phase diagram of the bulk material 

Before discussing the properties of the heterostructure, 
we will first analyze the magnetic phase diagram of the 



0.15 




0.00^ 



0.00 0.25 



0.50 0.75 1.00 
n 



FIG. 2: Phase diagram of the two-orbital model on the 4x4x4 
cubic lattice with J|f ^ = jf ^ = J^^ and A = 1.5. The resuhs 
are obtained comparing the energies of the G-AFM, C-AFM, 
CE-AFM, E-AFM, A-AFM, and FM states. 



corresponding bulk material. On one hand, this allows 
us to determine the model parameters to be used for 
the calculation of heterostructures; on the other hand, a 
better knowledge of the bulk phase diagram also helps 
in understanding the possible spin structures in the het- 
erostructure. 

Previous theoretical investigations22i have shown that 
in the two-orbital model the bandwidth depends on the 
electron-phonon coupling strength A and the superex- 
change coupling strength J^^. The larger the A and 
J^^ are, the narrower the bandwidth is. In Fig. [5] the 
phase diagram of the two-orbital model for the bulk is 
shown at various electronic densities, using A = 1.5 and 
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J . The phase diagram is obtained by 
comparing energies of several candidate states: A-AFM, 
G-AFM, C-AFM, E-AFM, CE-AFM, and FM states on 
the 4x4x4 lattice. Note that there could exist even more 
exotic states in the phase diagram, such as the CxE,i-x 
state previously proposed^^ at n > 0.5 if larger lattice 
sizes could be considered, and spiral states at n = 1 when 
in the presence of spin frustrationj^i, However, as shown 
below our interest will be mainly in the A-AFM state 
stabilized at electronic density n — 1, the CE and A- 
AFM states at n = 0.5, and the G-AFM state at n = 0, 
namely in regions where the C^T^i-x and spiral states are 
not expected to be relevant. Hence in the current study, 
the Cx^i^x and spirals states are not considered. 

The bulk phase diagram is very rich. At small su- 
perexchange couplings the ground state is FM in a broad 
density regime, as expected from the double-exchange 
mechanism. At larger superexchange coupling the sys- 
tem transitions, from low to high electronic densities, 
from G-AFM, to C-AFM, to CE-AFM, and finally to E- 
AFM phases, respectively. Close to n = 0.5 the CE phase 
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FIG. 3: Layer-averaged electronic density n{Z) and electro- 
static potential 4>{Z), at A = 1.5 and J^^ = 0.065, vs. layer 
index. Ua and nj, refer to the electronic densities of the a and 
b orbitals, and ntot = Ua + rib, tiq is the Fourier transform of 
the local electronic density in each layer at Q = (7r,7r). (a) 
are results at a = 1.0; (b) are results at q = 0.3. 



is stable over a wide range of J^^ values and has an al- 
ternate charge/orbital order; whereas the G-AFM phase 
at low n has neither charge nor orbital order. There 
are also two A-AFM phases: close to electronic density 
n = 1 at intermediate J^^ and close to n = at small 
J^^ values. The one close to n = 1 has the correct alter- 
nate Sx^ — r^/3j/^ — orbital order expected at n = 1 
from experimental information.— Thus, it is interest- 
ing to point out that in the range 0.05 ^ J^^ ^ 0.075 
the phase diagram consists of G-AFM, C-AFM, CE, FM, 
and A-AFM phases consecutively from low to high elec- 
tronic densities, correctly resembling the phase diagram 
of real narrow-to-intermediate bandwidth bulk mangan- 
ites,26 Hence, jf^ = jf^ = 0.065 and A = 1.5 will be 

used as couplings for calculations of the heterostructures 
to be discussed later in this section. 



B. CE state close to the interface of the 
heterostructure 

1. Emergence of CE properties near the interface 

In this subsection, the physical properties of the het- 
erostructure are investigated using the above described 
model parameters. In Fig. [31 the averaged electronic den- 
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FIG. 4: Optimized spin configurations of the heterostructure 
studied here, using A = 1.5, J*^ = 0.065, and a = 1.0. 
The dashed lines highlight the spin zigzag chains that are 
characteristic of CE states, with FM order within each chain. 



sity n{Z) and electrostatic potential (/)(.^) in each layer 
are presented at a = 1.0 and a = 0.3. At the two end 
layers of the heterostructure, the densities n{Z = 1) w 1 
and n{Z = 8) « converge to the expected values in 
the corresponding bulk materials. But charges are re- 
distributed in the rest of the layers due to the long- 
range Coulomb interactions. Note that there exists a 
plateau at ritot(^) ~ 0.5 in the layers close to the inter- 
face. To understand this feature, let us study the Fourier 
transform of the local electronic density in each layer, 
riQ = riie^^'^^ , where iVxy is the number of sites 

in each layer, nq at Q = (tt, tt) displays a peak in the 
layers close to the interface, where the plateau in ritot 
exists. This suggests the presence of a charge ordered 
phase with ntot(-Z') ~ 0.5 that is stabilized at the inter- 
face. Such a charge ordered state is found to appear at 
the interface for all the a values considered in our study. 
Hence, we believe the existence of this state is an intrinsic 
property at the interface of the heterostructure, at least 
within the approximations used in our calculations. 

To better understand the states that are located at or 
close to the interface, it is important to study how the 
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FIG. 5: Layer dependence of the spin structure factor S(k) for 
the optimized spin configurations shown in Fig. |4] at several 
momenta k: (0,0), (7r,0), (37r/2, 7r/2), and (7r,7r). 



spins are ordered in each layer of the heterostructure. We 
have observed that the spin order is not much sensitive to 
the value of a.^^ Thus, here only the information for the 
optimized real-space spin patterns at a — 1.0 is presented 
in Fig.m From this figure, it is clear that the layers close 
to the interface (layers 3, 4, and 5) exhibit a CE-type 
spin order. It is well known that in bulk systems and 
narrow bandwidth manganites, the CE state is stabilized 
at electronic density n « 0.5 and has a staggered charge 
order. Moreover, our previous results in this section using 
a cubic cluster to mimic the bulk also unveiled a CE state 
at the same density. Thus, the CE spin order found here 
further confirms that a charge/spin/orbital ordered state 
is stabilized at the interface.™ 

Moving from the interface towards one end of the het- 
erostructure (layer 1), the CE state gives way to the 
FM order in each layer. But spins in two adjacent lay- 
ers tend to be AFM coupled, see layers 1 and 2 for in- 
stance. Hence, the state at this end of the heterostruc- 
ture has an A-AFM tendency, resembling the spin or- 
der in the bulk. Similarly, the state on the other side 
of the heterostructure turns from the CE state to a G- 
AFM. The above described features are further confirmed 
by the layer-dependent spin structure factor S'(k) = 
Iv^ J2i j S» ■ Sj e*'' ('"*~'"j) shown in Fig. [5l Considering 
the degeneracy of states with ordering vector k = (tt, 0) 
and (0,7r), and with k = (7r/2,37r/2) and (37r/2, 7r/2), 
the CE order in layers 4 and 5 is nearly perfect. 



2. Novel states at the interface with no bulk analog in 
experimentally known phase diagrams 

The previous analysis shows that the interface has CE 
characteristics. This may be considered as an "obvious" 



result, since in a heterostructure the interpolation be- 
tween bulk materials with n = 1 and n = likely will 
induce n — 0.5 at the interface. In this simplistic con- 
ceptual framework, the properties at the interface can be 
guessed, with good accuracy, merely from the bulk phase 
diagram. While this provides a reasonable starting point 
to analyze results and make predictions, further analysis 
actually suggests that this is not the end of the story, 
and some surprises can be unveiled at interfaces. 

To illustrate this point, consider for instance layer 3. 
Here, the average density is very close to 1, yet the spins 
form zigzag chains as in the CE state of n = 0.5. The 
reason is that the spins at layer 3 are already being in- 
fiuenced by the robust CE state formed at layer 4. Thus, 
layer 3 is an exotic interpolation between the extremes 
case of the FM layers of the A-AFM state in one end, and 
the CE state at the interface. Moreover, note that the 
relative orientation between the spins of adjacent zigzag 
chains of layer 3 is not antiferromagnetic, as in a nor- 
mal CE state, but it has some canting. This "canted 
CE state" is not present in bulk phase diagram, to our 
knowledge. Note that a state called the "pseudo CE" was 
previously discussed in experimentsi^ Considering just 
the individual layers, this state is the same as the usual 
CE, but the coupling between CE layers is ferromagnetic, 
instead of antiferromagnetic as in the standard CE state. 
Thus, this pseudo CE state is not the same observed in 
our heterostructures. Also note that a "canted CE" state 
similar to that described in the present investigations 
was reported in early neutron diffraction experiments^ 
of Pro.rCao.sMnOa. However, further investigations for 
the same material22i2i suggested instead a mixed-phase 
interpretation of the results, with a mixture of AF and 
FM phases, as opposed to a uniform canted CE state. 
Thus, to our knowledge there is no evidence that the 
canted CE state exists in bulk form in real experiments, 
although more work is needed to fully address this mat- 
ter. Also we are not aware of previous theoretical inves- 
tigations reporting such a canted CE statci^ 

We also noticed that although the layers located right 
at the interface, i.e. layers 4 and 5, exhibit almost a 
perfect CE order individually, the spins in the two lay- 
ers are not perfectly AFM aligned with respect to one 
another compared with what they should be in a bulk 
CE phase. In other words, the CE state in the bulk 
is well known for showing a "stacking" property along 
the direction perpendicular to the CE plane, which is re- 
spected here with regards to charge and orbital but not 
with respect to the AFM relative order of the spins. In 
fact, in the layers Z=A and 5 of the heterostructure, the 
relative spin orientation is at approximately 90° degrees, 
namely they are spin perpendicular to one another. Once 
again, to our knowledge such an arrangement does not 
exist in the bulk. While it is obvious that the interplane 
spin deviation from AFM order in a CE arrangement will 
increase the superexchange energy of the system, these 
spin arrangements also allow for charge transferring out- 
of-plane. In other words, having an AFM order in a 
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Z=3 



Z=4 



Z=5 



Z=6 



FIG. 6: Real-space orbital pattern of the eg electrons for the 
layers close to the interface. The model parameters are the 
same as in Fig. ID The radian part of the electron wave func- 
tion (shown) is proportional to the local eg electronic den- 
sity m. The plot displays a transition from the staggered 
— /3y^ — order in layer 3, to a CE-type orbital order 
in layers 4 and 5, and then to a weak 3z^ — r'^ order in layer 
6, consistent with the spin patterns in Fig. U 



link reduces the effective hopping amplitude to zero in 
a double-exchange context. But if the spin order is not 
AFM, the kinetic energy in that link improves. Then, the 
system can gain double exchange energy in the z direc- 
tion at the expense of superexchange if the CE stacking 
does not involve AFM order. As a combined effect, the 
total energy appears to be lowered by this mechanism at 
the interfaces we studied here, while in the bulk it does 
not occur. 

The charge distribution in Fig. [3] suggests that on 
average the a {x^ — y"^) orbital has a higher occupa- 
tion number than the h (3z^ — r^) orbital. To ob- 
tain the exact orbital pattern at each site, the expec- 
tation values of the local pseudo-spin operators (rf ) and 
(rf) are calculated. Defining an effective phase angle 

= TT + tan^-^((Tf )/(r/)), we introduce a dressed state 

\b) = -sin(^j/2)c|a + cos(^j/2)4b |0> from which the 

orbital occupation is computed as {b\ni\h). The details 
of this type of calculations are well-known and they can 
be found in Ref. IH. 

To discuss the orbital pattern explicitly, let us focus 
on the optimized configuration with model parameters 
A = 1.5, J^^ = 0.065, and a = 1.0. For other a values 
the patterns look very similar. The orbital patterns in 
layers 1 and 2, where the bulk- like A- AFM phase is sta- 
bilized, display a clear staggered — r^/Sy^ — order. 
Reciprocally, layers 7 and 8 do not show any orbital or- 
der due to the vanishing value of the electronic density. 
However, the orbital patterns close to the interface in the 
range of layers from 3 to 6 are complicated and they are 
presented in Fig. [5] Here, we observe a transition from 
the staggered 2>x^ — jZy^ — order to the orbital order 
of the CE state, and then to a (very weak) 3z^ — or- 
der in the G-AFM state, with increasing layer index. In 
layer 3, the orbital pattern is very close to that expected 
of a n = 1 state, but note that the orbital population 



along one of the orientations of the diagonals is not iden- 
tical for each diagonal. This is caused by the influence 
of the CE state of the layer 4. In addition, in the CE 
state of layers 4 and 5, the "bridge" sites, which have 
a higher electronic density than the rest, have staggered 
Zx^ — jiy^ — order, but the "corner" sites show a 
uniform x^ — order, which is larger than in the CE 
state stabilized in the bulk. In general, we find that in 
layers where either FM or CE spin order exist the eg 
electrons prefer to form an in-plane orbital order. This 
is because in these layers the two oxygens connected to a 
Mn ion tend to shrink along the z direction to minimize 
the energy of electron-phonon interactions by decreasing 
Q3. But in the layers where there is a G-AFM state, the 
oxygens will expand along the z direction to partially 
compensate the shrinking effect along this direction in 



other layers.^ Hence, a 3z 



order may appear. In 



summary, the features of the orbital arrangements are 
dominated by what we expect to find in n — 1 A-AFM, 
n = 0.5 CE, and n = 0.0 G-AFM states. However, subtle 
deviations can be observed: non equivalent diagonals in 
layer 3, corner population in the zigzags of layers 4 and 
5, and weak orbital occupation along the z-axis in layer 
6. 



IV. RESULTS FOR WIDE BANDWIDTH 
MANGANITES 

In the previous section, the physical properties of the 
LaMnOa/CaMnOs and PrMnOs/CaMnOa heterostruc- 
tures were discussed. States with CE characteristics were 
found to be stabilized near the interface. But for the 
cases of the LaMnOa/SrMnOs and PrMnOa/SrMnOa 
heterostructures, we expect to obtain different interfa- 
cial states because the corresponding bulk materials, such 
as (La, Sr)Mn03 and (Pr, Sr)Mn03, have a wider band- 
width and the CE state is not stabilized in these com- 
pounds at half doping. 



A. The bulk phase diagram 

As in the previous section, let us first study the phase 
diagram of the wide bandwidth manganite in bulk form, 
focusing here on the case A ~ 1.2. This coupling is 
smaller than the one used for narrow bandwidth man- 
ganites. However, for J^^ — J^^ — J^^, as shown in 
Fig. [Tlja), the A-AFM phase at A = 1.2 cannot be sta- 
bilized near n « 1. This is because the A-AFM phase 
only appears in a very narrow regime of the phase di- 
agram in the two-orbital model if a cubic lattice sym- 
metry is considered. ^■^ This is compatible with the fact 
that the real compounds exhibiting the A-AFM phase, 
such as LaMnOa, have an orthorhombic, instead of cu- 
bic, lattice symmetry, with the lattice constant along 
the c axis smaller than those along the a and b axes<^ 
To better incorporate the lattice distortion found in the 
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FIG. 7: Phase diagram of the two-orbital model on the 4 x 
4x4 lattice with A = 1.2. Shown are the cases: (a) J^^ = 
jAF ^ jAF^ jAF ^ 2jfP/3 = J^^. The results 

were obtained comparing the energies of the states G-AFM, 
C-AFM, CE-AFM, A-AFM, E-AFM, and FM. 




FIG. 8: Layer-averaged electron density n{Z) and electro- 
static potential 4'{Z), obtained using A = 1.2 and J"^^ = 0.07. 
Ua and Uf, refer to the electron densities of the a and b orbitals, 
and ntot = na + nb. (a) are results for a = 1.0; (b) are results 
for a = 0.2. 



symmetry due to the strain effect. Thus, we adopt the 
inter-layer superexchange coupling = J^^ = 0.07 at 
the interface. 



bulk parent compound, an inter-layer superexchange cou- 
pling larger than the intra-layer one should be used, i.e. 
jfF > J|fF.ii^ The ratio Jf^/Jf^ is estimatedH>2i^ 
to be 1.2 '-^ 1.5. Using these numbers, the phase dia- 
gram for ^ 2jf '^/a = J^^ at A = 1.2 is shown in 
Fig.ElJb). The A-AFM phase is now stabilized in a wider 
regime of the phase diagram, including n k 1. Thus, 
tuning J^^ to 0.07, the system experiences transitions 
involving the G-C-A-FM-A phases with increasing elec- 
tronic density n, and this resembles properly the experi- 
mentally observed phases of wide-bandwidth manganites 
(R,A)M0.26 

As for the RMO/AMO heterostructure, the state on 
each side, far from the interface, must converge to its 
bulk phase. To obtain a stable A-AFM phase in the 
bulk RMO, we use the couplings A = 1.2 and J^^ — 

2jf F/3 J^^ = 0.07 on the RMO side of the het- 
erostructure. On the AMO side, since the bulk AMO 
still has the cubic lattice symmetryfi^. we adopt A = 1.2 
and J|^F = = J^^ = 0.07. At the interface, the 
lattice constants along the a and h axes will compress 
but those along the c axis elongate to recover the cubic 



B. The state at the interface of the heterostructure 

The averaged electrostatic potential and electronic 
densities in each layer ai a — 1.0 and a — 0.2 are pre- 
sented in Fig. [5] These charge distributions are actu- 
ally similar to those shown in Fig. [3) Interestingly, there 
is also a plateau at n « 0.5 indicating the existence of 
a fairly stable half-doped state near the interface. Al- 
though there is no CE phase in the phase diagram of 
wide bandwidth manganites, we find instead that the A- 
AFM phase can be stabilized at n « 0.5 in the bulk limit. 
Thus, the plateau at n fa 0.5 in Fig. [8] suggests the state 
at the interface to be the A-AFM state. This assumption 
is fully supported by the optimized real-space spin con- 
figuration results presented in Figs. [H] and [TOl for both 
a = 1.0 and a = 0.2, the A-AFM phase spin arrangement 
is found in layers close to the interface. 

Figures [5] and [TO] indicate that the transition from the 
A-AFM state at the interface to the G-AFM state at one 
end of the heterostructure depends on the value of the 
screening parameter a. At a large a value, an interesting 
result is found. In this case, there is one layer (layer 6) 
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FIG. 9: The optimized real-space spin configuration at each 
layer of the heterostructure with the model parameters used 
in Fig.iga). 



FIG. 10: The optimized real-space spin configuration at each 
layer of the heterostructure with the model parameters used 
in Fig.lHb). 



with an intermediate exotic state, that has no analog in 
the bulk experimental phase diagrams, to our knowledge. 
This state consists of alternate FM and AFM stripes, in- 
dicating a local mixed phase tendency. Note that we have 
carried out numerical studies on square clusters, A = 1.2, 
and J^^ = 0.07, simulating bulk two dimensional sys- 
tems, and in this case we do find a similar mixed AF-FM 
state at electronic quarter-filling densities. Thus, it is 
conceivable that this state may exist in bulk single-layer 
manganites at large hole doping as well. At small a val- 
ues, on the other hand, the transition is via introducing 
spin canting in several layers. We note that in real het- 
erostructures, the mixed-phase tendencies and spin cant- 
ing may coexist. 

Let us consider now the orbital occupation near the in- 
terface. In Fig. [Tl] the orbital pattern near the interface 
for A = 1.2, J^^ = 0.07, and a = 1.0 is presented. It 
is interesting to observe that the orbital pattern shows a 
clear transition from a staggered — r^/3y^ — order 
at one end of the heterostructure to a uniform — 
order at the interface, layers 4 and 5, although the spin 
order is unchanged and fixed into an A- AFM state. This 



is understandable, since we know in the bulk the A- AFM 
exists both close to n = 1 and close to n — 0.5. Com- 
paring with Fig. [Hl^a) , we observe that the orbital order 
is tightly connected to the average electronic density of 
the layer: the 2>x^ — /2>y^ — order appears at n « 1, 
but the — order is present at n ~ 0.5. Such a uni- 
form — y^ orbital order is also observed in the A- AFM 
phase at n « 0.5 in the bulk materiali^^iS Since the bulk 
A-AFM phase at n « 0.5 is metallic, the A-AFM state at 
the interface of the heterostructure can be anticipated to 
be a two-dimensional metal. However, note that such a 
2D metallic state could be unstable due to Anderson lo- 
calization introduced by the roughness and defects at the 
interface. Hence, an insulating behavior is more likely to 
be observed in real materials. 

In the exotic layer 6 exhibiting the mixed AFM-FM 
tendency, the orbital pattern is complicated. Accord- 
ing to Fig. [21 there are majority spins (black arrows in 
layer 6) and minority spins (red arrows). The minor- 
ity spin sites correspond to 3z^ — order. The sites 
adjacent to the minority spin sites display alternating 
Bx^ — r^/Sy^ — order, but the sites diagonal to the 
minority spin sites show a — order. Once again. 
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Z=3 Z=4 Z=5 Z=6 



FIG. 11: The orbital pattern found in layers 3 to 6 of the 
heterostructure. A uniform — orbital order is observed 
at the interface. The orbital order in layer 6 is exotic, i.e. not 
found in the bulk, similarly as the spin arrangement of the 
same layer. 



we remark that such an exotic orbital order has not been 
observed experimentally in bulk materials to our knowl- 
edge, although it may be part of theoretical phase dia- 
grams of models for two dimensional manganites in the 
bulk at large hole densities, and it is also conceivable that 
real single-layer manganites in the bulk may present also 
a similar phase. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, the magnetic and electronic properties of 
the states near the interface of RMO/AMO heterostruc- 
tures were investigated using numerical optimization 
techniques on small clusters at zero temperature. The 
states stabilized at the interface are found to be similar 
to the states in the bulk compound at electronic density 
n « 0.5. This is easy to understand. Let us consider 
the superlattice (RMO)j^/(AMO)„. When m,n > a/ a 
where a is the lattice constant along the z direction, there 
must exist layers exhibiting properties of the bulk parent 
compounds RMO and AMO. In this case, the electronic 
density in the regime near the interface is not much sen- 
sitive to m or n, but is always close to 0.5 if a is about 
the same on both sides of the heterostructure. Then, the 
regime close to the interface (within a/a in the z direc- 
tion) can accommodate the state stabilized at n « 0.5 
in the bulk. For the case of a wide bandwidth mangan- 
ite, a uniformly — ordered A-AFM phase is stabi- 
lized at the interface. Hence we propose that this state 
should exist in the long-period (LMO)in/ (SMO)n super- 
lattice where m and n are large enough such that there 
are bulk-like A-AFM and G-AFM insulating regimes far 
from the interface. It is quite reassuring that recent ex- 
periments on the (LM0)n/(SM0)2n superlattice at n = 3 
provided strong evidence for the existence of this orbital 
ordered A-AFM phaseJ^ However, for short-period su- 
perlattices, the electronic density near the interface may 
deviate from 0.5, and also note that a well-defined bulk- 
like regime may not exist in these structures. Hence the 
state stabilized at the interface may be different from 



the one in the long-period superlattice. For instance, 
in the (LaMn03)2n/(SrMn03)j^ superlattice with n < 4, 
the electronic density near the interface is always higher 
than 0.5fi^ falling into the FM regime in the phase di- 
agram of bulk LSMO. Hence, it is natural to observe a 
metallic FM state stabilized at the interface i^iSiiiiii^ii^ 
Increasing the number of Sr layers, bulk-like insulating 
regimes appear and the superlattice is driven through the 
MIT to be an insulator. The electronic density at the in- 
terface is also reduced by increasing the number of Sr 
layers. As shown in Ref. [l4|, the state at the interface is 
still FM with a 3z^ — orbital order. However, further 
increasing the number of Sr layers, the electronic density 
at the interface will approach 0.5. Then, the state at the 
interface becomes an A-AFM with a uniform x'^ — or- 
bital order as discussed above, and the properties of the 
heterostructure are dominated by bulk-like regimes.— 

When studying the wide bandwidth manganites, we 
have used J^^ > ■ For consistency, the same su- 
perexchange coupling ratio should be also used for the 
calculation of narrow-to-intermediate bandwidth man- 
ganites. But the A-AFM phase at n = 1 is already sta- 
bilized in a cubic lattice for a large electron-phonon cou- 
pling. Thus, changing the superexchange coupling ratio 
does not modify the phase diagram crucially. Calcula- 
tions have shown that setting J^^ > J^^ also does not 
change the results for the heterostructures. Hence, we 
have only presented here the results with J^^ = J^^ in 
Sec. mil In this case, interfaces dominated by CE-AFM 
characteristics have been observed. 

It is very important to remark that several of our re- 
sults have unveiled phases that are not observed in the 
bulk phase diagrams. They correspond to interesting 
modifications of the well-established bulk phases. For 
instance, for narrow bandwidth manganite heterostruc- 
tures, the existence of exotic spin arrangements, such 
as "canted CE" and others, have been reported in our 
investigations. For wide bandwidth manganites, unex- 
pected mixtures of FM and AF features were also iden- 
tified. While our observations obtained on small clusters 
need to be confirmed by other many-body techniques, 
the present computational studies revealed the possibil- 
ity of finding new phases at interfaces, that do not exist 
in the bulk. This is an exciting result that deserves fur- 
ther investigations. Figure [T^ summarizes schematically 
our results. 

Regarding size effects in our simulations, certainly it 
is possible that the actual spin, charge, and orbital pat- 
terns of the novel states may be more complicated than 
found in our present study that was limited to small sys- 
tems. However, note that in the heterostructures dis- 
cussed here the n = 1 A-AFM on one side, the n = 0.5 
CE or A-AFM at the interface, and the n = G-AFM on 
the other side are very robust and likely will be present 
in real heterostructures. With these states (with differ- 
ent layer electronic density) anchored somewhere in the 
structure, then it is very reasonable to expect magnetic 
states interpolating between them, that will have novel 
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FIG. 12: Schematic representation of the resuhs found in our 
investigations. I and V correspond to regions with properties 
similar to those of the bulk of the two materials involved in the 
heterostructures. Region III is very close to the interface. It is 
in this regime that the electronic density is approximately 0.5 
if the bulk components have densities 1 and 0, as in our study. 
Depending on bandwidths here either a CE or an A-AFM 
state are found. Finally, in regions II and IV, the material 
must interpolate in properties between I and III or III and 
V. These interpolations lead to states that do not appear to 
have bulk analogs. 



densities, then the novel states discussed here could also 
appear at the interfaces between puddles of the compet- 
ing phases as well. 

In conclusion, we find that the properties of the 
RMO/AMO heterostructure are closely, but not entirely, 
associated with the phase diagram of the bulk compound 
(R, A)MO. We summarize our main results in Fig.[T2l As 
one sees, although a uniform state does not exist in the 
heterostructure due to the redistribution of the charges, 
as a first approximation the state near the interface can 
be "read" from the phase diagram of the bulk compound 
(R, A)MO at electronic density n « 0.5, i.e., a CE state if 
(R, A)MO is a narrow-to-intermediate bandwidth man- 
ganite, but an A-AFM state with uniform — orbital 
order if (R, A)MO is a wide bandwidth manganite. How- 
ever, the intermediate states in between the interface and 
the bulk-like regimes are sensitive to both the model pa- 
rameters and the length of the heterostructure. At least 
in our studies, it can be either a spin canted CE state or a 
state showing local FM-AF mixed tendencies. We believe 
our results can be used to describe the ground state prop- 
erties of the long-period {YUAO)^/ [KMO)^ superlattices 
as well. 



properties. In the bulk when the electronic concentration 
is the same in every layer of course it never happens that 
a magnetic layer must interpolate between others with 
different densities. Thus, we are confident that novel 
magnetic states, of a form likely even more complex than 
unveiled here, would be present in manganite superstruc- 
tures if simulations using larger clusters were possible. 
Note that for bulk phase-separated manganites,™ involv- 
ing a competition between states with different electronic 
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